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The Bose distribution for a gas of nonrclativistic free bosons is derived in the 
framework of qp-deformed second quantization. Some thermodynamical functions 
for such a system in D dimensions arc derived. Bose-Einstein condensation is 
discussed in terms of the parameters q and p as well as a parameter v' which 
characterizes the representation space of the oscillator algebra. 

1 Introduction 

The use of quantum groups and quantum algebras is now largely displayed in 
theoretical physics. From an elementary point of view, quantum algebras may 
be thought of as q- (or qp)-dcformations of Lie algebras. Roughly speaking, 
there are two types of applications of g-deformations to Physics (cf. Ref. 1). 
Applications of type A dot not rely on phenomenology. For example, applica- 
tions of type A concern the quantum inverse scattering method, the quantum 
Yang-Baxter equation and a true definition of space-time. On the other hand, 
applications of type B are entirely phenomenological. Along this vein, we may 
mention the use of the two-parameter quantum algebra U qp (xi2) (a Hopf al- 
gebra) to rotational collective dynamics of atomic nuclei. 2 For applications of 
type B, we immediately face the problems of the physical interpretation and 
the nonuniversality of the deformation parameter(s). 

The considerations above can be applied to deformed oscillator algebras. 
Indeed, oscillator algebras and quantum algebras may be connected since the q- 
(or c^-deformed creation and annihilation operators, which are defined in an 
oscillator algebra, may serve for constructing realizations of q- (or <7p)-deformcd 
Lie algebras. 

It is the aim of this lecture to examine the interest of using qp-deformcd 
creation and annihilation operators in statistical mechanics of bosons in ar- 
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bitrary dimension. In other words, we ask the following question: What are 
the implications, at the level of the Bose distribution and the Bose-Einstein 
condensation, of replacing ordinary boson operators by gp-deformed creation 
and annihilation operators? Furthermore, is it possible to test a deformed Bose 
distribution and to find a physical interpretation of the deformation parame- 
ter^)? 

This subject was already investigated, mainly in the case of q-deformed 
boson operators, by many authors. 3-18 The case of gp-deformed boson oper- 
ators was considered in Refs. 17 and 18. Most of the works in Refs. 3 to 18 
use the Fock representation 19-25 for the deformed oscillator algebra inherent 
to the considered creation and annihilation operators. In contradistinction, 
the representation introduced by Rideau 26 (see also Ref. 1) for the g-deformed 
oscillator algebra is used in Ref. 16. 

We shall deal in this paper with both the Fock and the Rideau represen- 
tations of deformed oscillator algebras. These representations are described in 
Sec. 2. Then, Sec. 3 is devoted to gp-deformations of the Bose distribution 
and its consequence for the Bosc-Einstcin condensation phenomenon. Finally, 
some conclusions are presented in Sec. 4. Two appendices close this paper. 
In particular, Appendix B concerns a gp-deformed correlation factor for the 
radiation field. 

This work takes its origin in Refs. 9 and 18. It is however to be noticed 
that some of the developments and the conclusions given here substantially 
differ from the ones in Refs. 9 and 18. 

2 Deformed Oscillator Algebra 

Following Refs. 1 and 26, we characterize the algebra W qp q by three operators 
a, a + and N such that 

• the operator a + is the adjoint of a, 

• the operator N is self-adjoint with a discrete nondegenerate spectrum, 

and 

• the operators a, a + and N satisfy the commutation relations 

[N,a]=-a, [A,a+]=a + , aa+ - Qa+a = — !— [q N (q - Q) - p N (p - Q)} 

q-p 

(1) 

where the three parameters q, p and Q are fixed parameters taken (a priori) 
in the field of complex numbers. 

The algebra W qp Q spanned by the deformed annihilation (a), creation (a + ) 
and number (N) operators is called a deformed oscillator algebra. In the special 
cases Q — q and Q — p, the operators a and a + are referred to as qp-boson 
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operators. They entail the particular cases of the g-boson operators as used 
in mathematics (when p = 1) and in physics (when p = q^ 1 ). Four particular 
cases are of interest for the applications (see also Refs. 19 to 26). 

1) The case Q = q, p = 1: This case corresponds to 

[N, a] = -a, [N, a + ] = a + , aa + - qa + a = 1 (2) 

i.e. to g-bosons as used in the mathematical literature (see also Refs. 19 to 21 
and 25). 

2) The case Q = q, p = q^ 1 : This case corresponds to 

[N,a] = -a, [N,a+]=a + , aa+ -qa + a = q- N (3) 

i.e. to g-bosons as used in the physical literature. 22 ' 23 

3) The case Q = q^ 1 , p = q^ 1 : This case corresponds to 

[N, a] = -a, [N, a+] = a + , aa+ - q^a+a = q N (4) 

i.e. to g-bosons as used in the physical literature. 22 ' 23 

4) The case Q = 0: This case corresponds to 

N+l _ N+l 

[N, a] = -a, [N, a + ] = a + , aa+ = + £ (5) 

q-p 

which admits two interesting subcases corresponding to either p = 1 or p = g _1 . 

The nondeformed oscillator algebra corresponds to (Q — q, p = 1, q — > 1) 
or (Q = q, p = q — > 1) or (Q = q -1 , p = q — > 1). Note that the cases 
(Q = 0, p = 1, q — ► 1) and (Q = 0, p = q -1 , q — ► 1) correspond to the exotic 
limiting case where aa + = 1. 

The Hilbertean irreducible representations of the oscillator algebra W q = 
W qq -i q were dealt in details by Rideau. 26 The Hilbertean irreducible represen- 
tations of the more general algebra W qp Q can be derived equally well. It is 
enough to say here that the relations 

a\n) - ( Q"A + p"° °" ~ P " + q"° g " ~ °" V \n - 1) 

/ — r>™ +1 n n + 1 — D n + 1 \ i 

a+\n) - (q^X 0+p ^ ^ + q»° q - *«—) \n + l) 

N\n) = {v Q +n)\n) (6) 

provide us with a formal representation of W qp Q . We thus obtain an infinity of 
Hilbertean irreducible representations of W qp Q that may be arranged into two 



types: R{v'q) with n € N and S(Xo, Vq) with n £ Z. We shall deal in this work 
only with the (Fock) representation R(0) of the algebra W qpq and the (Fock) 
representation R{y' ), with v* Q ^ 0, of the algebra W q . We shall not use the 
case of the (nonFock) representation S(Xq, vq). This representation which has 
no limit a,t q = p = Q = 1 is briefly discussed in Appendix A. 

The (deformed) Fock representation R(0) of W qpq is obtained from Eq. (||) 
by taking Xq = v§ = and n G N. Then, we have 



a\n) = y[[n]} qp \n-l), 
a + \n) = \J[[n+l]] qp \n + l), 

N \n) = n \n) (7) 
with a |0) = 0. In Eq. (0), we use the notation 

M\ qP = (8) 

q — p 

It is to be observed that the operator N is a nondeformed operator that co- 
incides with the usual number operator (in the sense that its spectrum is N). 
Note that 



} g ,P=i = ^— j, [[x]] 9J>=q -i = g g _ pi CM 



The case of the nondeformed (i.e. usual) Fock representation is reached for 
(p = 1, (jf — s- 1) or (p = q — > 1). Finally, the case (p = 1, q = 0) 

corresponds to the exotic situation where 

a\n) = |n-l), a + \n) = \n + l), N \n) = n\n) (10) 

Going back to the representation i?(0) (given by Eq. (]?])) of W qpq , we observe 
that 

a+a = [[N]) qp , aa+ = [[N + l)] qp (11) 

so that the relation 

aa + -pa+a = q N (12) 
holds for this representation in addition to the defining relation 

aa + — qa + a = p N (13) 

Equations ( |l2| ) and (^3|) clearly show that the representation R(0) of W qpq 
exhibits the q «-> p symmetry. 



The (deformed) Fock representation R[v'q) of W q is obtained from Eq. 
by taking Ao = 0, vq = v' ^ and n£N. Then, we have 



a+ \n) = q '*y/[n+l] q \n + 1) 
N\n) = K + n)|n) (14) 
with a |0) = 0. In Eq. (|l4]), we use the notation 

Mq = [[x]]q,p=q-i (15) 

and we have q > 1 and v' € R. It is to be observed that the operator N does 
not coincide with the usual number operator; the usual number operator is 
now N — v'q (in the sense that its spectrum is N). As an important result, we 
have that if g — ^ 1 + then v' Q — > 0; in other words, we can write 

lim R{v' ) = usual Fock representation (16) 

q^l + 

For the representation R(v' ), we have 

a+a = q-<[N - v' ] q , aa + = q~ v '° [N + 1 — v' ] q (17) 

Consequently, we verify that the relations 

aa+ - qa+a = q~ N , aa+ - q^a+a = q N ~ 2u '« (18) 

hold in W q . Equation ( |Ts| ) indicates that the representation R(is' ) of W q does 
not present the q <-> q~ l symmetry. 



3 Bose Distribution and Bose-Einstein Condensation 

Let us consider a gas of nonrelativistic free bosons. Its Hamiltonian (in a 
second-quantized form) reads 

H:=Y,H k , H k := (E k - y)v k (19) 

fc 

In Eq. @, \x is the chemical potential while E k and v k are the kinetic energy 
of a boson and the number operator for the bosons, in the k mode, respectively. 
The Bose factor for the k mode is then 

fk ■■= |tr(e-" H a+ 0jb ) (20) 
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where 

Z:=tr(e- pa ) (21) 

is the partition function and (3 = (fcsT) -1 the reciprocal temperature. 

It is essential to note that Eqs.(|l9|)-(pT|) have the same form as in the 
nondeformed case. They correspond to the basic formulas for the statistical 
mechanics of bosons. Here, we do not deform the basic formulas. The trace and 
the exponential in ( p0[ ) and ( pl| ) are nondeformed functions. The deformation 
is introduced via the use of annihilation (eifc) and creation (a£) operators that 



satisfy gp-dcformed commutation relations. The form chosen for Hk in Eq. (19) 
is nothing but the common one where the usual number operator v% takes its 
eigenvalues in N. The deformation cannot enter the theory at this level by 



replacing z/& by a^a k in (19) since the operator a^a k has not its eigenvalues in 
N. On the contrary, we must keep a^a k rather than i>k in (p0[). Indeed, should 
we use Vk in place of a£a k in (ppf), we would get a nondeformed value of the 
Bose factor. As a conclusion, Eqs. (fL9l)-(f2l|) constitute a reasonable starting 
point for a (minimal) deformation of the Bose distribution. 

We are now in a position to derive several types of deformed Bose distri- 
butions. The deformation is explicitly introduced by using either the represen- 
tation i?(0) of W q pq (see Eq. (R)) or the representation R(v' ), with v' Q ^ 0, of 
W q (see Eq. (|l4|)). The two representations can be treated (in a formal way) 
on the same footing by considering (cf. Eqs. ([?]) and (14)) the representation 



flfe I"-} = P 2 \/[[n}] qp \n - 1) 
+ 



X, \ n ) = P 2 V [[«+!]]«> 



Nk\n) - K + n)\n) (22) 

of the algebra W qpq . We assume (mathematical hyphothesis) : {dfe,a£ ,Nk} 
and {aj,a^,Nj^ are two commuting sets when j ^ k. Furthermore, we take 
(physical hyphothesis): 

v k := N k - u' Q (23) 



As a trivial result, we have 



Z 



where 



V := /3{E k - n) (25) 
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The partition function Z is thus independent of the deformation parameters q 
and p. Then, the Bose factor fk can be written 

f k = p v 'o{l - e-")tT(e-^ N ^o)[[ Nk _ v ' Q ]] qp) (26) 

which can be shown to converge in each of the following five cases: 

• (q,p) e R+ x R+ with < q < e"^ and < p < e - ^ 

• (q,p) 6 R + x R+ with p = q^ 1 and c 13 ^ < q < e~ /3/i 

• (q,p) e R + x R+ with p = 1 and < q < e"^ 

• (q,p) e C x C with p = q and < \q\ < e _/3A1 

• {q,p) £ S 1 x S 1 with p = q = q~ x . 

The latter two cases, although acceptable from the mathematical point of 
view, do not lead to physically acceptable statistical distributions. 
As a formal result, Eq. ( Pq) yields 

\q — p ei — q p — q e 11 — p J 

or alternatively 

In the limiting situation where q — ► 1 and p — > 1, we recover the classical Bose 
distribution 

* fe = e^T (29) 
The distribution can be rewritten as 



fk = P v ° 



1 $fe , p- 1 $ A 



9 -p 1 + (1 - q)<S> k p-q\ + (l-p)$ k 



(30) 



in term of $>k. 

For the purpose of practical applications, it is useful to know the develop- 
ment of the distribution /& in integer series. In this respect, we have 

h = p< V e-"^ +1 ) ( q —qi + *^±pi) (31) 

p5 \q-p p-q J 

which reduces to the expansion Y^jLi e ~ v ^ °f the ordinary Bose factor when 
q — * 1 and p — ► 1. 

In the thermodynamical limit, the energy spectrum of the system of bosons 
may be considered as a continuum. Thus, fk is replaced by the qp-dependent 
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factor /(e) which is fk with Ek = s. Therefore, in physical applications, we 
have to consider integrals of the type 



poo 

J s := / e s f(e)de 
Jo 



(32) 



By using the development (plf), Eq. (32) leads to 



J s =T{s + l)(k B T) s+1 a{s + l) 



Ql" 



S > -1 



(33) 



where 



/3m(j'+i) 



a(s + l) qp = p< (j + 1)s+ i ([[j + !]]«- 



j=0 



igp) 



(34) 



In Eq. (|33]), T is the Euler integral of the second type. In the case where /x = 0, 
we have a(s + l) qp — > C( s + 1) f° r <7 ~^ 1 an d p — > 1, where £(s + 1) is the 
Riemann Zeta series (that converges for s > 0). 

We now derive some thermodynamical quantities, in D dimensions, for a 
free gas of N(D) bosons contained in a volume V(D). The density 



p{D) :- 



N(D) 
V{D) 



(35) 



is given by 9 



p(D) = N (D)T(-)(k B Tpa(-) 



(36) 



where 



iV (£>) 



D 



m 2 



2(2tt)t r(f + 1)" n 1 



(37) 



In Eq. (|37|), 5 is the degree of (spin) degeneracy and m the mass of a boson. 
The total energy can then be calculated to be 



E(D) = N (D)V(D)T(^ + l)(fc fl T)T+V(^ + l) 



(38) 



The specific heat at constant volume easily follows from Cy = (f^)y We 
obtain 



C V (D) = ^N(D)k B 



2 Jqp 



d 



(- 

l2+ j 2 -(#-l) 



(39) 



Furthermore, the entropy is 



S(D) = N(D)k B 



+ K 



\ 2 / qp 



knT 



Note that the state equation for the gas of bosons is 

2 E(D) 



p(D) 



DV{D) 



(40) 



(41) 



so that the pressure p(D) assumes the same form as in the nondeformed case. 

Finally, other thermodynamical quantities may be determined in a simple 
manner by using the thermodynamic potential 



n(D) 



D 



N (D)V(D)J R 



D 



E(D) 



(42) 



where Jn is given by ( p3| ) and 

We now examine the condensation of a system of bosons in D dimensions. 
The corresponding density for such a system is given by (|3^) which can be 
rewritten as 

p(D) = N Q (D)J^_ 1 (43) 
as a function of the integral Jo_ _, . 

2 1 

As in the classical case (i.e. q = p = 1), we have to define a procedure for 
generating the critical Bose temperature T C (D). For this purpose, we follow 
Ref. 9 by taking 

(3fi = j3/zi — lng, q > 1 (44) 

where [X\ is a negative constant corresponding to the chemical potential in the 
classical case. This g-dependence of the chemical potential /i makes it possible 
to have a good behaviour of the distribution f% when /ii — » _ only for the two 
situations p = q^ 1 and p — 1. Therefore, in the remaining part of this section, 
the parameter p stands for either p = q~ or p = 1. The Bose-Einstein conden- 
sation is then obtained by introducing /xi = in Eq. d43|). Thus, we find that 
the critical temperature below which we obtain Bose-Einstein condensation is 
given by 



1 



T c (D) qp = p-o»o 



P(D) 



1 



N (D) T(§)a (§) qp 



where 

"<T> 



^ (i + i)" 



for 



P = 9 



and 1 < g < e 



-/3m 



(45) 



(46) 
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or 

D 1 00 1 
M^)<,p := - V -w for p = 1 and 1< g < e"^ (47) 

2 ? J=0 + !) 2 

We thus end up with three models: a two-parameter model Mi with (v' 6 R*, 
p = q^ 1 , 1 < q < e - ^) and two one-parameter models, viz. M2 with (v' = 0, 
p = q^ 1 , 1 < q < c^ 13 ^) and M 3 with (u' = 0, p = 1, 1 < g < c - ^). For each 
of these models, the Bose-Einstein condensation phenomenon takes place only 
when the dimension D is greater than 2, exactly as in the classical case. (The 
series in Eqs. (||) and Q) converge for D > 3.) The models Mi, M 2 and M 3 
can be easily tested for D = 3 by comparing for each model the theoretical 
and experimental temperatures T c (3) for 4 He superfluid in phase II. For the 
models M2 and M3, we obtain 

T c (3) qp > T c (3)classical > 7c(3)cxp ~ 2.17 K (48) 

Therefore, the models M2 and M3 give results which are worse than the one 
given by the classical model (corresponding to q = p = 1) for which we have 



Tc (3)classical — ^c(*^)(/=l,p=l 



mkB 



P(3) 
2.612 3 



(49) 



Unlike the models M 2 and M 3 , for the model Mi, it is possible to find couples 
(v'o,q) or which T c (3) q p=q -i is in agreement with the experimental value. As 
a point of fact, our test rules out the models M 2 and M3. 



4 Conclusions 

Among the various gp-deformations of the Bose factor studied in this work, 
only the ones corresponding to p = q^ 1 or p = 1 yield acceptable statistical 
distributions for bosons. However, there is no (/-deformation, with u' — 0, of 
the Bose factor that leads to a correct value of the Bose-Einstein transition tem- 
perature. This has been tested for 4 He but we note that T c (3) qp > T , c (3) c iassicai 
holds in general. The latter result generally applies to other q-deformations, 
available in the literature, 4,5 ' 10 ' 15 ' 16 based on the use of the representation R(0) 
of the oscillator algebra W q . 

The model Mi deserves more optimistic concluding remarks. As a matter 
of fact, for this model, corresponding to (v' < 0, p — q~ x , 1 < q < e _,3M ), it is 
possible to decrease the critical temperature T c (D) qp and to get Bose-Einstein 
condensation temperatures in accordance with experiment. However, as a 
drawback, the model Mi depends on two parameters v' and q. The parameter 
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q brings a translation factor (3 1 \nq for the chemical potential since Eq. (44) 
gives 

M = Mi = Mi - k B Tlnq (50) 

that clearly shows that the importance of the g-deformation increases with the 
temperature. The parameter v' is more difficult to interpret. This param- 
eter is essential, via the factor qn v a in Eq. (jig), for decreasing the critical 
temperature. In this respect, it might be of importance for describing some 
interactions between bosons. 

To close this paper, it is worth noting that alternative choices, non-linear 
in iVfc, are possible for the Hamiltonian H in Eq. (|l9|). In this direction, we 
may mention, in the case p — q^ 1 , the works of Refs. 4, 5, 10, 15, and 16. 
We hope to return on these matters in the future, especially on the difficult 
problem of the interpretation of the parameter v' . 
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Appendix A 

We deal in this appendix with the representation S(Xq,Uq) of the oscillator 
algebra W q with < q < 1. This representation follows from 

a\n) = (<f Ao + ^N,)' |n-l) 

a+\n) = ( 9 " +1 Ao + + \n+l) 

N |n) = (v + n) \n) (51) 

where n £ Z and the real parameters q, Xq and vq satisfy the constraint 
equation 

g —^_ <Ao<5 -o _L_ (52) 
1-q 1-q 

A characteristic of this representation is that the spectrum of the operator N 
is not bounded from below. It is thus hardly feasible to connect iV to the usual 
number operator. Another important characteristic concerns the absence of a 
limit of ^(Aoj^o) when q — > 1~. In fact, the generic representation S(\q,uo) 
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disappears at q = 1 since Eq. ( ]5^) becomes meaningless for g — > 1~. The 
fact that representation S(\q, vq) has no limit at q = 1 is probably the reason 
why this representation has never been used (to the best of our knowledge) in 
physical applications. 



Appendix B 

We discuss here the consequence of a gp-deformation of the correlation function 
g( 2 ) of order two associated to the radiation field. In the nondeformed case, 
we know that g( 2 ' takes two values, viz. g^ = 1 for a coherent monomode 
radiation and g^ = 2 for a chaotic monomode radiation. An interesting 
question arises: Is it possible to interpolate between the latter two values by 
replacing the ordinary boson operators by gp-deformed boson operators? In 
this appendix, we limit ourselves to the representation i?(0) of the oscillator 
algebra W qpq . 

Our basic hypothesis is to describe the radiation field by an assembly of 
photons = 0) with the Hamiltonian 

h := hk, hk := huJkVk (53) 
k 

The corresponding Bose statistical distribution for the k mode is then given 
by ( f28[) where 77 is replaced by £ := /3/kJfe- In the case of gp-bosons, we adopt 
the definition 

g m := (54) 

where a and a + stand for the annihilation and creation operators for the 
k mode. The abbreviation (X) in ([54]) denotes the mean statistical value 
Z~ l ti (e~ l3h X s ) for an operator X. Equation ( |5^ ) can be developed as 

(55) 

Of course, (a + a) = /j. as given by ( ]28| ) with r) = £. In addition, the other 
average values in (p5|) can be calculated to be 

l( + \2\ (e^ - l)(et + qp) 

« a a ) ) - {e ,_ q 2 ){e( _ qp){e ,_ p2) ( 56 ) 



and 



(q N a + a) = q— 6 ] (57) 

(e« - g 2 )(e« - qp) 
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Finally, we obtain 



g ^ = {q+p) ^L l w J^^g z _ (58) 

with convergence conditions that parallel the ones for (with — 2/3 replaced 
by £). [For instance, when G R + x R + , we must have < q < and 

< p < ea.] The qp-deformed factor g&> depends on the parameters q and p 
in a symmetrical manner (<?<-> p symmetry). It also presents a dependence on 
the energy huik of the k mode and on the temperature T. In the limiting case 
where q — > 1 and p — > 1, we get (j*- 2 -* = 2 that turns out to be the value for the 
choatic monomode radiation of the black body. Let us now examine the cases 
of low temperatures and high energies. In these cases, e^ is the dominating 
term in each of the differences occuring in Eq. (^8|). Therefore, we have 

g^-q + p (59) 



at low temperature or high energy. Equation (|5S|) shows that we can interpolate 
in a continuous way from g^ 1 = 1 (coherent phase) to — 2 (chaotic phase). 
It is thus possible to reach the value g^ — 1 without employing coherent states. 
It should be observed that we can even obtain either g^ > 2 or g^ < 1 from 
Eq. (|58|). The situation where g^ < 1 may be interesting for describing 
antibunching effects of the light field. 
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